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Abstract 

The unintegrated gluon density in the photon is determined, using the Kimber-Martin- 
Ryskin prescription . In addition, a model of the unintegrated gluon is proposed, based on 
the saturation model extended to the large-x region. These gluon densities are applied to 
obtain cross sections for charm and bottom production in and 77 collisions using the kt 
factorization approach. We investigate both direct and resolved photon contributions and 
make comparison with the results from the collinear approach and the experimental data. 
An enhancement of the cross section due to inclusion of non-zero transverse momenta 
of the gluons is found. The charm production cross section is consistent with the data. 
The data exceed our conservative estimate for bottom production in collisions, but 
theoretical uncertainties are too large to claim a significant inconsistency. A substantial 
discrepancy between theory and the experiment is found for 77 bbX, not being cured 
by the kt factorization approach. 



1 Introduction 



Production of heavy quarks at high energies has been vigorously studied experimentally in 
recent years. Measurements of cross sections for charm and bottom production have been 
performed in proton-proton proton-photon and photon-photon [0,^] collisions. 

The charm data in 7*p and 77 collisions may be reasonably well described by the standard 
collinear formalism, based on LO QCD with NLO corrections 0. For bottom however, the 
experimental results exceed significantly the theoretical expectations in all cases. The largest 
discrepancy has been found for bottom production in 77 collisions at LEP p[, where the 
measured cross section is larger by a factor of four than the QCD prediction. 

The enhancement of the bottom production cross section was reported with different beams 
and at different energies, which suggests the presence of an important systematic effect, omitted 
in the QCD analysis. It is particularly puzzling because of the large mass of the bottom quark, 
giving a safe ground for the perturbative approach. Attempts have been made |T0|, pA|, p!2|, p!^. 



IJ, |T5|, |T^, [T^ to resolve this problem by going beyond the standard collinear formalism 
and use the kt factorization approach |]1O|,0,[T^,|T9[. Thus, instead of assuming that massless 
partons are distributed in the colliding objects having a negligibly small transverse momentum, 
one considers the complete kinematics of parton scattering, taking into account the transverse 
momenta. 

This intrinsic transverse momentum k of the parton is built up in the perturbative evolution, 
as a result of subsequent emissions of gluons or quarks and its distribution is parameterized 
by the unintegrated parton distribution. The influence of the parton transverse momentum 
on the cross section depends on the relevant hard matrix element, which has to be evaluated 
for virtual partons (off-shell matrix element). Calculations using the off-shell matrix elements 
combined with the unintegrated parton distributions were performed for bottom production in 
pp collisions |l^,|T3|,|15| . Indeed, the obtained cross sections are larger than those in the collinear 
approximation and agree with the data within uncertainties. For 7*p an enhancement is also 
found for a direct photon |T6|,P, 



TBIl, but is not sufficient to describe the data. 



The unintegrated gluon distribution in the proton evaluated at the factorization scale 
J-'g{x, k^, /i^) is a subject of intensive studies itself (for a review, see |^). This quantity depends 
on more degrees of freedom than the collinear parton density, and is therefore less constrained 
by the experimental data. Various approaches to model the unintegrated gluon have been 
proposed. For instance, in the leading logarithmic l/x approximation, evolution of J-'g{x, k , /z^) 
is given by the BFKL |^T| or CCFM equations. The unintegrated gluons following from 
those equations were fitted successfully to inclusive data from ep scattering ||2^, |4| . This 
approach is restricted to the small x regime. Recently, it has been shown p5|,^^ that the 
information contained in the collinear parton densities combined with the properties of parton 
emission amplitudes (e.g. the angular ordering) is sufficient to determine the unintegrated gluon 
up to a large x. An interesting model for the gluon is also given by the successful saturation 



model [^, introduced by Golec-Biernat and Wiisthoff (GBW). 

Models for the unintegrated parton distributions in the photon were not available until very 
recently [^,^^ and no results for the resolved photon are known beyond the collinear limit. 
Thus, the main purpose of this study is to obtain in an independent way, the unintegrated 
gluon distributions in the photon, using the Kimber-Martin-Ryskin (KMR) prescription 
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and to apply them to describe the heavy quark production (charm and bottom) in 77 and jp 
coUisions. Apphcation of the kt factorization formahsm for the case of resolved photon(s) is 
performed for the first time. We also obtained an alternative gluon density in the photon based 
on the generalized saturation model ||32| for 77 interactions. We will explore a variety of gluon 
parameterizations and account the other model ambiguities in order to estimate the theoretical 
uncertainties. We will examine whether the excess of the bottom production in these processes 
can be explained within the kf factorization approach. 

The paper is organized as follows: in Section |^ the unintegrated gluon in the photon is 
obtained and its properties are discussed, in Section ^ the kt factorization formulae are presented 
and in Section ^ the cross sections for heavy quark production in 77 and 7p collisions are 
calculated. A discussion of the results is given in Section ||, followed by conclusions in Section ||. 



2 Unintegrated gluon distributions in the photon 

In the construction of unintegrated gluon distributions in the photon we apply the same method 
as for the distributions in the proton. The off-shell parton distributions in the proton are better 
known and their properties have been investigated to a great extent (see ||2^ and references 
there- in), while similar distributions in the photon are poorly known and no attempts have 
been made to describe them until recently p9|,pO|,|3T||. 



2.1 The KMR approach 



The conventional gluon distribution g{x, jj?) corresponds to the density of gluons in the photon 
having a longitudinal momentum fraction x at the factorization scale /x. This distribution 
satisfies the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution [33] in /i^ and it is 
universal for the photon in different processes. The distribution does not contain information 
about the transverse momenta k of the gluon, which is integrated over up to the factorization 
scale 11 



xg[x, /i 



(X 



(1) 



However, in order to better describe processes by properly considering the transverse momentum 
of the gluon, unintegrated parton distributions J^g{x,k'^ , jj?) are needed. These distributions 
take into account the complete kinematics of the partons entering the process at the leading 
order (LO). 

As a first attempt, the unintegrated gluon density may simply be obtained [25], at very low 
x, from the coUinear gluon density xg{x,^'^). This density becomes independent of the hard 
scale /i^, and will only depend on one scale fc^ 



J- g[x, k , /i 



d 



dQ- 



■[xg{x,Q'^)] 



r=k' 



(2) 



The above equation no longer holds, as x increases, since Tgix, k^.Q"^) becomes negative. 
This may be circumvented, however, by introducing a Sudakov form factor Tg{Q,fi), which 
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takes into account subleading corrections at low x. Thus, the unintegrated distribution has 
now a 2-scale dependence [^,|27[| 



d 



[xg{x,Q^)^Tg{Q,^i)] 



(3) 



with the form of Tg{Q, /i) given below. The form factor represents the probability of the gluon 
with the transverse momentum k to survive untouched in the evolution up to the factorization 
scale. 

A better framework for unifying the small x and large x regions is provided by the Catani- 
Ciafaloni-Fiorani-Marchesini (CCFM) equation |^2|]. This equation is an evolution equation 
for the unintegrated gluon distribution J-'g(x, , fi^) which considers real gluon emission in a 
ladder and is based on angular ordering of the gluons in the chain. The formalism has a natural 
interplay of two scales: the transverse momentum k of the gluon and the hard scale /z, which 
corresponds to the maximal angle of emitted gluons. Thus, the unintegrated gluon distributions 
which can be constructed will have a 2-scale dependence J-'g{x, fc^, fi^), where n will have a dual 
role, that of factorization scale and controlling the angular ordering. At small x, the CCFM 
formalism is equivalent, in the leading log(l/a;) approximation, to the Balitskij-Fadin-Kuraev- 
Lipatov (BFKL) formalism [21], and J-'g{x, which satisfies the BFKL equation becomes 



/^-independent. At moderate x, the angular ordering is replaced by fcjordering, and the CCFM 
equation reduces to DGLAP. 

A simplified solution to the complicated 2-scale CCFM evolution was obtained by Kimber, 
Martin and Ryskin (KMR) in [^]. They observed that the /idependence in the distributions 
enters only in the last step of the evolution, and single-scale evolution equations can be used 
up to the last step. In this approximation, the unintegrated gluon distribution is given by 



J- g{x, k , /i 



Tg{k,fi) asik^ 



k' 



2n 
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where Pgg{z) and Pgq{z) are the gluon and the quark splitting functions, while g{x,Q^) and 
q{x, Q"^) are the conventional gluon and quark densities. The Sudakov form factor introduces 
the dependence on the second scale /xin the last step of the evolution and has the following 
form 



Tg{k,fi) = exp 



k' 



2-n 



1-5 



dzz 



qgy 



(5) 



where 6 = -p^ is chosen to provide the correct angular ordering of the real gluon emissions. 

We have extended the KMR formalism to the case of the photon. In the following we will 
use the unintegrated gluon density in the photon defined by equation (^. The conventional 
gluon and quark densities in the photon are expressed following the Gliick-Reya-Schienbein 
(GRS) parameterization Thus, g{x, Q'^) and g(a;, Q'^) consist of two components, a pointhke 
(perturbative) component parameterized in ||35| and a hadronic component, given by the parton 
distribution functions in the pion I^Uj. For instance, for the gluon we use 



xg{x,Q'^) 



1 



-xgpi{x, Q2) + (^2 + Gl)xg^{x, Q^), 



(6) 
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a) b) 

Figure 1: Comparison between the off-shell gluon distributions in the proton and in the photon, 
using KMR approach, as a function of (a) longitudinal momentum fraction x and (b) transverse 
momenta k of the gluon, for fixed values of (GeV^) and x, respectively. 



with = 0.043 and = 0.50, while the respective formulae for quarks can be found in |35 
The obtained unintegrated gluon density is defined for k"^ > = 0.5 GeV^, which is the 
starting scale for the GRS distribution. However, an extrapolation to cover the whole range in 
k^ has been performed, extending the gluon density to values of fc^ < k^hy normalizing it to 
the GRS distribution in the following way J-'g{x, k'^^jj?) = xg{x, kD/k^. 



Another solution to the CCFM equation was found using the "single loop" approxi- 
mation, when small-x effects can be neglected in the CCFM equation for medium and large 
X. Thus an exact analytic solution can be obtained, and a comparison between this analytic 



solution for the proton and the KMR approximation shows quite good agreement Simi- 
larly, for the photon the unintegrated gluon distributions obtained from the exact solution 
of the CCFM equation in the single loop approximation can be well represented by the KMR 
distributions constructed using the integrated quark and gluon distributions and the Sudakov 
form factor. 

Although the KMR constructions of unintegrated gluon distributions for the photon and 
proton are similar, the distribution in the photon is significantly different due to the pointlike 
component. As can be seen in Fig. |l], a direct comparison between the unintegrated gluon 
in the photon and in the proton shows a relative enhancement for large x and large k in 
the case of the photon. This enhancement is due to gluon emissions from the perturbative 
quark box, making the gluon distribution much harder as compared to the proton for large 
values of x. For small values of x, the similar shape of both distributions indicates that the 
information about the shapes, contained in the input at large x, is partially lost in evolution. 
Note that for the respective KMR gluon density in the proton we have used the conventional 
GRV parameterization [RM] for the proton. 
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2.2 The GBW gluon 



Another parameterization of the unintegrated gluon density in the photon can be obtained 
using a simple generalization of the Golec-Biernat and Wiisthoff (GBW) parameterization of 
the gluon in the proton |2^. The unintegrated gluon density introduced in for the proton 



:F,{x,e) = ^Rl{x)k'exp{-Rl{x)e), (7) 

is related to the effective dipole cross section within the saturation model 

a{x,r') = ao[l -exp(-rV(4i?^(x)))], (8) 

which describes the interaction between a proton and a color qq dipole coming from a photon 
fluctuation 

Nf ^ 

a^*p = y [ dz I (fr \^>\z,r)\^a{x,r^). (9) 

In the above equations, r denotes the transverse separation of the quarks and z gives the 
longitudinal momentum of the quark in the photon. The wave function of the photon is 
represented by r), where a indexes the flavor of the quark in the dipole, and its form can 
be found in [^, for instance. Thus, for real photons the mass of quark a gives the characteristic 
scale for the dipole distribution in the photon. The saturation radius is given by 

with ao = 29.12 mb, a, = 0.2, Qo = 1 GeV, A = 0.277 and xq = 0.41 x 10"^. The three free 
parameters ctq, A and xq have been fitted and the saturation model describes successfully both 



inclusive and diffractive 7*p scattering [^]. The mass of the light quarks u, d and s, m, 



0.21 GeV is taken from a fit of the saturation model to inclusive two-photon observables p2[ . 

Following the generalization of the GBW saturation model for the case of 7*7* scattering 
introduced in [^], one can easily construct in a similar fashion the unintegrated gluon distri- 
bution in the photon. In such a case, we consider the scattering of two color dipoles, into which 
the photons fluctuate, one light qq dipole and one heavy QQ dipole 



1 1 
dz, I rfVi|vl/"(;2i,ri)p/^ dz2 I rfV2|^«(z2,r2)paf(x,ri,r2). (11) 



The heavier dipole with the transverse separation r2 provides the hard scale at which the 
dipole content of the second photon is probed. In this configuration, the relative size of the 
heavy dipole ((^2) ~ l/2mQ) is smaller than that of the light dipole ((^i) ~ l/2mq). The 
effective dipole-dipole cross section a'^'^{x,ri,r2) is a generalization of the GBW cross section 
from eq. (|), introducing an effective dipole separation radius r^ff depending on the size of the 
two dipoles [^]. In our configuration with one heavy dipole r| ^ rj, the effective cross section 
reduces to eq. (Rl), depending only on r2. 
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Figure 2: The gluon distribution xg{x) in the photon as a function of x. The data points 
describe the gluon content in the photon extracted from the photoproduction of hard dijets 
(meanp^ = 74 GeV^) observed at HERA [^. The GRS parameterization |^ is compared with 
the integrated GBW distribution from (|l^ at fi^ = Pt, with and without the phenomenological 
threshold factor included. 



In this region, the integrals from the equation (|rT|) over Zi and Vi of the first dipole can be 
performed independently. In the leading logarithm log^mQ/nig) approximation, the result of 
this integration is dominated by the contribution 



Nci{fi) = J2 (fr^\^\z^,r^)\\ (12) 

with a lower cut-off in Vi provided by the typical size of the heavy dipole, ~ l/2mQ. This 
integral may be interpreted as the number of dipoles in the photon at the scale given by the 
mass of the heavy quark. The final result after the integration will be a form for the 7*7* cross 
section which is similar to the 7*p cross section in eq. (|]) 



Nd{2mQ) X f dz2 I ciV2|^^(^2,r2)|' a(x,r2). (13) 



The number of dipoles Nd in the photon is found to be 1.46-aem for charm production (Mc = 
1.3GeV) and 2.43-Q;em for bottom production (Mf, = 4.8 GeV). 

The extraction of the gluon density in the photon from eq. (|13D will give a parameterization 
similar to the one in eq. (^ 

^,(x,fc^/i2) = N,{f,) X ^^Rl{x)eexp{-Rl{x)e). (14) 

This includes a multiplicative factor given by the number of dipoles in the photon, and the 
parameter do = |(Toas introduced in the generalization of the GBW model (the factor 2/3 
is a reminder of the quark counting rule, with ctq representing the cross section in the blackness 
limit for the photon, and respectively, do for the proton). The number of dipoles that enters in 
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a) b) 

Figure 3: Comparison between the unintegrated gluon distributions in the photon, using the 
KMR approach (dashed hues) and the generahzed GBW parameterization (sohd hues) as a 
function of (a) longitudinal momentum fraction x and (b) transverse momenta k of the gluon, 
for fixed values of (GeV^) and x, respectively. 



the gluon density has an intrinsic dependence on the hard scale given by the heavy quark mass, 
which propagates as a secondary scale at the level of the unintegrated gluon distribution. All 
other parameters from the original GBW parameterization are kept unchanged. 

As described in in order to extend the color dipole model to moderate and large 



X values, the introduction of phenomenologically motivated threshold factors was necessary. 
Thus, we have imposed the following form on the total cross section in 77 interactions, a'^*'^ ~ 
(1 — In the kt factorization approach, this factor can be understood as having a 

dual contribution, from the off-shell matrix elements and the unintegrated gluon distributions, 
which consider the correct kinematics of the hard process. To account for the full kinematics 
in the unintegrated density of the photon, we will introduce such a multiplicative factor in this 
density. When probing the gluon content of a hadron with a photon, only sea quarks can be 
picked, so the number of spectators is Uspect = 4in the case of a proton (3 constituent quarks 
plus 1 sea quark), and Uspect = 3 for the photon case (2 quarks from the dipole plus 1 from the 
sea). 



Multiplying J-'g{x,k'^, fj?) from ([IID with the factor (1 — x)^ and integrating it to obtain 



the corresponding on-shell gluon density, we find the same dependence with x as the GRS 
distribution for large-x. As seen in Fig. 0, the integrated gluon distribution provides a better 



agreement with the existing data extracted from photoproduction of hard dijets at HERA ||39 
as compared to the case where no threshold factor was included. Similarly, we have extended 
the applicability of the GBW gluon distribution for the proton (|^) for large values of x by 
introducing the multiplicative factor (1 — xY . 

The above unintegrated gluon distribution obtained using the extended saturation model 
exhibits different x and fcdependence than the previous density stemming from the KMR ap- 
proach. These differences are best expressed in Fig. ^, where a suppression for large values of 
k and an enhancement at small momenta can be seen in the GBW gluon. In spite of their 
unlikeness, the two distributions give quite similar results when integrated over the transverse 
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momenta. Figure ^ shows how the integrated KMR distribution, similar to the conventional 
gluon density GRS, and the integrated GBW gluon compare with data. Let us note how poorly 
constraining data is for the gluon content of the photon, and that a new fit could give an in- 
crease in the gluon distribution which can alter significantly predictions for cross sections based 
on these distributions. 



3 Cross sections for heavy quark production 

Heavy quarks may be produced in two-photon collisions by one of the three mechanisms: a 
direct production (Fig. ^), a photopro duct ion off a resolved photon (Fig. ^b) and a by a double 
resolved process (Fig. Ijc). The direct contribution to the process 77 QQX is governed by 
elementary QED amplitudes |^D|. In the case of proton-photon scattering, this component is 
absent. 




a) b) c) 



Figure 4: Diagrams illustrating the heavy quark production in 77or 7p collisions through 
different mechanisms: a) direct production (only for 77 interactions), b) photoproduction in 
single resolved case, where one photon or the proton is resolved, and c) double resolved case, 
where both incoming particles are resolved. 



In the coUinear limit one uses the leading twist term of the operator product expansion, 
neglecting the transverse momenta of the partons. The cross section for heavy quark photo- 
production off one of the photons being resolved reads 

^tUs,M') =Y.I -/.(x,/x^)a^{(M^x.), (15) 
i •' ^ 

where the partonic cross sections d'^{ are well known up to the NLO approximation and 
fi{x , fi'^)is the parton i distribution function in the photon, at the factorization scale /i ~ 2M. 
In the above, s denotes the 77 invariant mass squared, and M is the quark mass. 
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In the kt factorization formalism, the complete kinematics of the gluon-photon fusion is 
taken into account and the small light-cone component of the longitudinal momentum of the 
gluon is integrated out. Then, the cross section takes the following form 

= E / -—H^.k\^^')a'4{k\M\xs\ (16) 

in which the unintegrated parton density k'^^fj?) and the off-shell partonic cross section 
{k^ ^ M'^ ^ X1X2S) are employed. The partonic cross sections are evaluated using off-shell 
matrix elements. The form of a^^ is well known in the literature, see for example |TT|, and we 
quote it in Appendix 0. Formulae (|15D and (|16D , with the appropriate substitution of parton 
densities, are valid also for the photoproduction off the proton. It is important to note that, 
in a 77 collision one or the other photon may be resolved, thus the inclusive cross section for 
heavy quark production acquires an additional factor of 2. At the LO approximation, only 
gluons contribute to the QQ production. 

The double-resolved contribution to the process 77 QQX is described by 

dXi dX2 ^ ,,2^^Y^_ ,,2^;^c/^ 

Xx X2 



^Us.M^) =J2[ — — /.(xl,;u2)/,(x2,/i2)(7^/(M^a;lX2.), (17) 

• • J X-\ Xo 



in the collinear limit and the cross section in the kt factorization formalism reads 

dxi dx2 d^ki d^k2 

Xi X2 



M') = J:[ ^^^^^ kl /i^) ^,(X2, kl ^^') at/ik,, k2, M\ x,X2s), (18) 

"j^ J Xi X2 IT TT 



where af/ for gluons is given in Appendix P (following from I^TJ). Analogously, one of the 



photons may be replaced by the proton in order to obtain a resolved photon contribution to 
heavy quark photoproduction off the proton. 

In the following, we shall restrict ourselves to the effects of transverse momentum in the 
gluon kinematics, and the quark contribution will only be taken in the collinear approximation. 
This should not affect significantly the results, as the heavy quark production is driven mostly 
by exchanges of gluons. For the collinear limit, all the theoretical estimates were obtained using 
the PYTHIA Monte Carlo p|. 



4 Results for heavy quark production 

The measurements of cross sections for the inclusive charm and bottom production in e^e~ 
and ep collisions were performed at LEP and HERA respectively. For bottom production at 
LEP, only the total rate e^e~ — >■ e^e~hhX was measured |p, representing an average of the 
77 — ^ cross section weighted with the flux of photons in the electrons. For charm, the cross 
section (7(77 — >■ ccX) was determined for different 77 collision energies |]^. At HERA, the cross 
section for 7p bbXis known ^ for the collision energy averaged between W = 94 GeV and 
W = 266 GeV. The data for charm production at HERA [|| are available for virtual photons, 
at different virtualities and collision energies W in the form of F2^°'^"^{x,Q'^). We give in 
this section the theoretical estimates for these processes based on the kt factorization formalism 
and study the theoretical uncertainties. 
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4.1 Theoretical uncertainties 



The cross sections for heavy quark production are described by formulae (|ly) and ( |T8|) with the 
LO partonic cross sections given by ( p3D and (pTI). The results of the numerical evaluation of 
these formulae depend somewhat on the model for unintegrated parton densities and the choice 
of parameters. We have examined the following options for different elements of the model: 



The heavy quark mass M is plagued by a fundamental uncertainty due to confinement of 
color - there are no free quarks, and consequently, there is no on-shell quark mass. The running 
quark mass in QCD varies with the scale. It is not clear, at which scale the quark masses in 
the matrix elements should be evaluated, because the virtualities of heavy quarks are different 
for different lines in the relevant Feynman diagrams. Thus, we shall assume for the h quark 
that 4.5 GeV < < 5 GeV and for the c quark that 1.3 GeV < < 1.5 GeV. 

The energy scale /2 that enters the running formula of the strong coupling constant 
in the partonic cross section (see Appendix) is usually chosen to be of order of the typical 
momentum transfer characterizing the process. However, the optimal value of /2 is such, that 
the contribution of higher orders in the perturbative expansion is minimal. Thus, without 
knowledge of higher order corrections, p, is uncertain and, in order to account for this ambiguity 
we considered three options: (1) Jj? = M'^+p^ (standard), (2) Ji^ = {M"^ +p^) /A (low scale) and 
(3) Jj? = 4M^ (large scale), where p = k (gluon transverse momentum) for the gluon-photon 
fusion (^) and p = ki — k2 for the two-gluon process (pT|). 

Running of the coupling constant. We use the standard one- loop running formula for ag, 
with four flavors. We default choice, Aqcd = 140 MeV, such that as(M|) = 0.117 

as given by the latest QCD fits. We also test the PYTHIA default value: Aqcd = 250 MeV, 
corresponding to the value of as(M|) = 0.128. 



Unintegrated parton distributions. For the proton, we take into account the CCFM gluon 
distribution, as given in the CASCADE Monte Carlo |H3[, the unintegrated gluon obtained 



from the GRV and MRST parameterizations using the KMR method (uGRV and uMRST) and 
the gluon following from the saturation model (GBW). In the case of a real photon, the KMR 
method is applied to the GRS gluon (uGRS) and an alternative model of the unintegrated gluon 
is given by the saturation model for photons (GBW7), as explained in Sec. |2.1| . Furthermore, we 
vary the factorization scale in the two-scale gluon distribution J-'g{x, fe^, yU^) between /i = Mand 



2M. 



Parton momentum fraction x. In the coUinear approximation, one neglects the non-zero 
transverse momentum k of the incoming parton whereas in the kt factorization approach, the 
transverse momentum is included in the kinematics of partonic scattering. Thus, for the virtual 
photon-gluon fusion, the invariant mass of the system is s^g = xs in the coUinear approximation 
and s'^g = xs — k when the complete kinematics are taken into account. In a standard approx- 
imation method, the cross section in the coUinear factorization can be obtained from the one 
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in the kt factorization, using a substitution a^^ {xs^k"^ .Q"^) a''f(xs,k'^ = 0,Q'^)9{Q'^ — fc^), 
and the integrated parton distributions as defined by the eq. (p. However, such approxima- 
tion neglects the fact that the kinematical threshold in DIS for a^-^ depends on fc^, i.e. the 
kinematical threshold for the virtual photoproduction is located at x ~ Q'^/s in the coUinear 
approximation and at 

in the kt factorization framework. 

In order to investigate the inclusion of the latter threshold treatment, the standard approx- 
imation could suffer modifications by introducing a rescaled variable z = xQ'^/ (Q^ + k"^)- This 
will improve the approximation and lead to an alternative relation between the kt factorization 
and collinear approximation expressions, where 

zgiz, = 1"^' TMQ' + k'')/Q\ k\ Q2). (20) 

J 

Thus, the use of the rescaled variable is a way to quantify the ambiguity in obtaining the 
unintegrated gluon distribution from the integrated one. 

In our investigation of the effects of the threshold treatment in the evaluation of the heavy 
quark production cross section, we consider a substitution 



1 + fcV(4M2) 

in the unitegrated gluon distributions jFg(a;', fc^, /i^). Note, that for massive quark photopro- 
duction, 4M^ replaces Q^. Such rescaling is not necessary for the CCFM gluon, where the 
kinematical effects of the transverse momentum are already accounted for in the F2 fits. 



(21) 



4.2 Results for interactions 

In Fig. ^ we give a set of results for cross sections for 7p — > hhX with the direct photon (Fig. |^a) 
and with the resolved photon (Fig. ^) . The default results (full curves) are obtained by taking 
the CCFM unintegrated gluon in the proton (from CASCADE), Mb = 4.5 GeV, /i = 10 GeV, 



Aqcd = 140 MeV and the standard scale for the running coupling (see Sec. |4.1| ). For the gluon 
in the photon, the KMR method is applied to obtain the unintegrated GRS parameterization 
(uGRS). Besides the default choice, we also use the unintegrated GRV distribution (uGRV), 
based on the GRV NLO parameterization, and the GBW gluonQ, leaving other parameters 
unchanged. Furthermore, we investigate the impact of the x-rescaling, for the uGRV gluon 
(direct photon) and for both the uGRV and uGRS gluon (resolved photon). We modify the 
default choice of the QCD parameters: we use Aqcd = 250 MeV and the low scale for running, 
in order to obtain an upper limit for the cross section. A very conservative estimate follows 
from the choice Mb = 5 GeV and the large energy scale for a^. We have checked, that using the 
unintegrated MRST gluon and variations of the renormalization scale n have a minor infiuence 
on the results, hence we do not include the corresponding curves in the figure. 

^In the definition of the unintegrated gluon in the saturation model, one assumes a fixed as = 0.2, and, 
consequently, the same choice was made for the coupling in matrix elements when GBW gluon was used. 
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Figure 5: Cross sections of bb production in 7p interactions from the kt factorization approach, 
showing a) the single resolved contribution and b) the double resolved contribution. Details of 
the curves are presented in the text. 



It is visible in the figure, that the contribution from the resolved photon is significant - 
roughly 20-30% of the direct photon cross section for the default choice of parameters. We 
have checked, that in all cases, the resolved photon contribution rises faster with the energy 
W than the direct photon one. The results are rather stable against modifications of the 
unintegrated gluon and the quark mass. The largest contribution to the uncertainty of the 
cross section comes from the details of the QCD running of as and the choice of energy scale 
- reflecting a potential influence of higher order corrections. When the low scale of a^and the 
large value of Aqcd are assumed, the direct contribution gets enhanced by about 50% and the 
resolved one gets doubled. In this extreme case, the calculation results for the sum of the direct 
and resolved photon cross sections a(W = 180 GeV)= 143 nb is close to the experimental data 
point a = 206 ± 19^40 (see Fig. |^). Thus, although the QCD predictions are significantly 
lower than the data, inconsistency cannot be claimed, due to the large experimental error and 
the large theoretical uncertainty coming from uncontrolled higher order corrections. 

The saturation model estimate of the direct photon contribution to bottom production is 
lower than the default result. The double resolved contribution obtained using the GBW pa- 
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Figure 6: Comparison between the collinear factorization and the /c^ factorization approaches to 
bb production in 7J9 interactions. Shown are the experimental data from EMC and HI ||5|. 
a) The total bb cross section in the kt factorization (solid line) includes single and double 
resolved contributions. In the collinear factorization, the total cross section (short dashed) also 
includes the specific contribution from the qq annihilation for the double resolved case (long 
dashed), b) Comparison between the single resolved components in the kt (solid line) and 
collinear (short dashed) factorizations, and between the double resolved components in the kt 
(dotted) and collinear (long dashed) factorizations. 



rameterization for both the photon and the proton is slightly larger than its KMR counterpart. 
The constraint we have imposed on the parton density at large x via the threshold factor plays 
a very important role for the resolved photon case, being dependent on the gluon at relatively 
large x. The consistency of the standard QCD and saturation model results is not surprising, 
as both sets of parameterizations are constrained by the same experimental data (see Fig. 
It is, however, striking that the unintegrated GBW gluon leads to a lower cross section than 
the uGRS gluon, while giving higher integrated distribution, as shown in Fig. 0. This effect 
is caused by strong suppression of the GBW gluon at transverse momenta larger than the 
saturation scale, as we demonstrated in Fig. |]b. 

In Fig. ^a we show results for the total cross section, a{'jp bbX) obtained in the kt 
factorization formalism and in the collinear approximation, compared to the experimental data. 
In both cases, the direct and resolved photon contributions are added. The resolved photon 
contribution, coming from the partonic process qq QQis only evaluated in the collinear limit 
and demonstrated to be negligibly small. In Fig. ^b, the cross sections are decomposed into 
the direct and resolved photon components. 

We used the unintegrated gluon in the proton obtained via the KMR method (uGRV) and 
the default values of the other parameters for the kt factorization calculation, as given in Sec. |4.2| . 
The collinear limit results are obtained using the corresponding integrated gluon distributions 
(GRV), transverse momenta are set to zero in the definition of energy scale for as, that is 
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Figure 7: The structure function F2'^ as a function of x at different values of between 
1.8 GeV^ and 130 GeV^. The continuous hues represent the results from the saturation model, 
with the parameters from |^|, compared with the experimental results from ZEUS Qj. 



we take as(M^). Other parameters take their default values. Therefore, we gain insight into 
the actual difference between the studied approximations, not caused by discrepancy of input 
parameters. Thus, it is clear from Fig. Pd, that a large enhancement (by a factor of three) of 
the resolved photon cross section occurs due to non-zero transverse momentum effects, whereas 
the direct photon contribution gains only modest 10-20% in the magnitude. The total cross 
section is larger by some 20-30% in the kf factorization approach. 

It has been checked that charm production at HERA is well described within the kt factor- 
ization approach [^. For completeness, we show in Fig. ^ our results for the charm structure 
function F2^"'^"^{x,Q'^), based on the standard saturation model, compared to experimental 
data. The quark mass Mc = 1.3 GeV. The agreement of the theory and the data is very good, 
in contrast to the discussed bottom production case. The results we have obtained for bottom 
production in the single resolved case are compatible with previous studies of heavy quark pro- 
duction using the kt factorization approach [l^,|T^. These studies found agreement with the 
first erroneous data point from HERA 0, which gave a cross section for bottom production 
lower by a factor 2 than found later ||^. With the inclusion of the resolved photon in the 
kt factorization framework, we restore a similar level of agreement with the experimental data. 



4.3 Results for 77 interactions 

The default set of parameters for bottom production in 77 collisions is the same as it was in 
the case of the 7p process. The unintegrated gluon in the photon is unfolded from the GRS 
parameterization, using the KMR method (see Sec. |2.1| ). We investigate the sensitivity of the 
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Figure 8: Cross sections of bb production in 77 interactions from the kt factorization approach, 
showing a) the single resolved contribution and b) the double resolved contribution. Details of 
the curves are presented in the text. 



cross section to variations of the parameters. The results are shown as a function of 77 energy 
W in Fig. p. The default results (continuous line) are compared with the results incorporating 
the kinematical rescaling of x (see Sec. Besides that, the low scale in the running formula 
of as and the large Aq^are assumed and the case of = 5.0 GeV and the large scale of 
as is shown. Furthermore, the unintegrated gluon in the photon from the saturation model 
is used. We stress that within the presented W range, the results are driven by the gluon 
at relatively large x ~ 0.1, where the saturation model is less reliable and phenomenological 
threshold factors need to be imposed. 

Furthermore in Fig. we illustrate the enhancement of the cross sections due to non- 
zero transverse momentum of the gluon. This figure is constructed in strict analogy to Fig. ^ 
discussed in Sec. [4.2| . The conclusions from these results are also very similar to those obtained 
in the jp case. Let us only mention, that the integrated and unintegrated GRS parameterization 
of the gluon in the photon was used. Note, that the QPM contribution (see the diagram in 
Fig. ^) shown in Fig. ^ is included in the total cross sections (7(77 bbX). In order to obtain 
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Figure 9: Comparison between the collinear factorization and the hi factorization approaches 
to hb production in 77interactions. a) The total 66 cross section is represented with the sohd 
hne in the kt factorization and with the short dashed hne in the colhnear factorization. The 
long dashed line gives the specific contribution from the qq annihilation in the double resolved 
case, corresponding to the collinear approach. The dotted line shows the contribution from the 
quark box (see Fig. ^a). b) Comparison between the single resolved components in the kt (solid 
line) and collinear (short dashed) factorizations, and between the double resolved components 
in the kt (dotted) and collinear (long dashed) factorizations. 



the total cross section a[e^e~ e'^e~bbX), the 77 cross section needs to be weighted with 
the photon flux in the electrons. Thus, we obtained the value of a{e~^e~ — e~^e~bbX) = 1.9 
pb for the default case (based on KMR gluon density), while the use of the GBW gluons 
gives a{e~^e~ e~^e~bbX) = 2.7 pb. The latter value is lower than what we previously 
found using a generalization of the saturation model [^], with the difference coming from the 
threshold factor imposed at the level of the unintegrated gluon ^^(a;, fc^) instead of the total 
cross section. The behavior of the cross section in the vicinity of the kinematic threshold 
is crucial for interpretation of the LEP measurements of bottom production. An interesting 
discussion of this problem may be found in |17]. 

The general picture, which emerges, may be summarized with the following. The theoretical 
uncertainty of the QCD results for bottom production is rather large — of the order of 50% for 
the single resolved case and, even larger for both photons resolved. The double resolved photon 
contribution is only a small correction (a few percent) to the single resolved photon over the 
studied energy range. The QCD models give results 2 times smaller than the GBW-model. 
The total cross section following from the kt factorization scheme are not significantly larger 
than the ones obtained in the collinear limit. Therefore, the enhancement due to the use of 
unintegrated gluon is not sufficient to solve the 6-excess puzzle in 77 collisions. 

In an analogous scheme as for the bottom, we present in Fig. ^ the results for (7(77 ccX). 
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Figure 10: Cross sections of cc production in 77 interactions from the kt factorization approach, 
showing a) the single resolved contribution and b) the double resolved contribution. Details of 
the curves are presented in the text. 



Because of the lower quark mass, the uncertainties related to the energy scale choice and 
the definition of gluonic x are larger than for the bottom production. The double resolved 
photon contribution is again small, few percent correction to the dominant single resolved 
photon contribution. The data points are fitted well within the uncertainty band of theoretical 
estimates, as seen in Fig. |Tl|. The extended GBW model as introduced in and QCD give 
similar results for charm within the shown energy range, whereas in the bottom case the GBW 
gluon gave a larger cross section than QCD. The main reason of this difference is that, within 
the considered range of energies, the gluon is probed at lower x for charm, as compared to the 
bottom case. 



5 Discussion 



The main goal of this study is to investigate whether the puzzle of bottom production excess 
in 7p and 77 collisions may be naturally explained in the kt factorization framework. Thus, we 
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have estimated the cross sections for direct and resolved photon(s) and compared the results 
to their counterparts in standard coUinear approximations. 

We, indeed, found some enhancement of the cross section for a direct photon scattering off 
both the proton and the resolved photon. This effect itself is too small, though, to get the theory- 
results close to the experimental data. Furthermore, we showed, that the cross sections driven 
by two-gluon fusion are substantially larger (by a factor of 2-3) when the unintegrated gluon 
is used. Despite this enhancement, these subprocesses contribute only as a 20-30% correction 
to the total cross section for bottom production in the 7p collisions, and less than 10% for 77. 

The sensitivity of the results to various model uncertainties is found to be large. Thus, 
in the marginal case, the theoretical results for 7p — »• bbX in the kt factorization framework, 
is enhanced and agrees with the HI data point within errors. We do not interpret this as 
a strong indication of consistency between data and the theory, but rather as a consequence 
of the wide uncertainty band. Better understanding of higher order corrections is crucial to 
determine whether the HI data point contradicts QCD results. 

The picture is much more clear for 77 — > bbX. This process is dominated by the QED 
box diagram at low energies and by the single resolved photon mechanism at larger energies. 
Irrespectively of the assumption made within our framework, the emerging results for e~^e~ — >■ 
e~^e~bbX, cr*^ = 1.9 pb at the e~^e~ collision energy ^/s = 200 GeV, are more than three 
standard deviations below the experimental data, 13.1 ± 2.0 (stat) ± 2.4 (syst) pb Thus, in 
the standard QCD+QED approach, the b production in e~^e~ collisions is hard to explain. 

Interestingly enough, the agreement between the theory and the experimental data is good 
for the charm case both for 77 and 7*p. The very different behavior of charm and bottom cross 
sections is surprising and calls for an explanation. 
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6 Conclusions 



In this study we have analyzed photoproduction of heavy flavors in 77 and 7p collisions using 
the kt factorization approach. First, we obtained parameterizations of unintegrated gluon in 
the photon using the KMR method. We compared features of the unintegrated gluon in the 
photon and in the proton. Furthermore, we proposed a parameterization of the gluon in the 
photon based on a generalization of the saturation model. 

The parameterizations, combined with off-shell matrix elements, were used to estimate cross 
sections for charm and bottom production, including contributions from resolved photon. The 
impact of non-zero transverse momentum of gluons was studied. Some enhancement of the 
cross section was found in the kt factorization approach. In particular, we demonstrated the 
importance of the resolved photon contribution to bottom production in 7p collisions. Sensitiv- 
ity of the theoretical estimates to the details of the model was investigated. The conclusion is, 
that the use of kt factorization approach brings the theoretical results for bottom production 
closer to the data, but large discrepancies remain. For 7p collisions, a major inconsistency 
cannot be claimed, because of large experimental errors and theoretical uncertainties. For the 
77 case, the 6-production excess is statistically significant, despite the uncertainties. The kt 
approach, based on QCD, does not agree with the data from LEP at the three sigma level. On 
the contrary, the charm production is well understood within QCD. 

Thus, an interesting question arises, why the 6-quark production excess is found in various 
processes. A potential explanation may be, perhaps, provided by higher order corrections. 
Still, it is not clear why similar, or even more important, corrections would not affect charm 
production in a similar way. It is also interesting to ask about uncertainties of parton densities 
in the photon since the experimental constraints from measurements of the photon structure 
and a jet photoproduction are not very stringent. Definitely, an attempt should be made to 
perform a new global fit of parton densities in the photon, including the bottom production data. 
Another interesting explanation of the 6-excess at the Tevatron was suggested in [^] where the 
non-perturbative fragmentation function of the b quark into B mesons was updated by fitting to 
precise LEP data. It was found, that the discrepancy between standard theoretical calculations 
and the Tevatron data is significantly reduced when the improved model of fragmentation is 
used. Finally, there is also a possibility, that the answer to the question leads beyond QCD |^ . 
Therefore, it is important to further constrain theoretical uncertainties of bottom production 
rate estimates in QCD. 
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A Virtual gluon-photon fusion 



The partonic cross section for an off-shell gluon with the transverse momentum fcdepend on 
the density matrix for gluon polarizations In this study, we assume that the gluons have 

the BFKL-like polarization tensor, that is 

E^^^r = ^. (22) 

The result of the integration over the phase for an off-shell matrix element describing a gluon- 
photon fusion 'y{p)g*{k) — >■ QQ is well known. We quote the result from [ pJ] 



with 



+ [8 + p - (2 + 3p)L(/?)]^ + [-8 + 2m] (^^^ I (23) 
i^ = 2p-k, p = , (24) 



\ 



l-p(l-^) ' (25) 



and 



The heavy quark charge is denoted by eg. 

B Cross sections for virtual gluons 



m = l\ogl±4- (26) 



In this appendix, we follow conventions of [^, with some minor modifications. The kinematics 
of the gluon fusion process g*{ki) + g*{k2) — > Qij>i) + Q{p^)i with respect to the four- vectors 
Pi and p2 of the incoming photons/hadrons, are given by 

ki = Xipi + ki, 

k2 = X2P2 - k2, 

Pz = (1 - zi)xipi + Z2X2P2 + fci - A, 

P4 = ZiXiPi + (1 - Z2)X2P2 - ^2 + A. (27) 



The two-body phase space of the QQ pair can be written in the following 



way 



,,<^..J^^?- /A.f.-^±Mi-,M. (28) 

87r^^i(l - zi) \ 2:1(1 -2:1) 
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where 



q = ki- k2, 
V = XiX2{2pi ■ P2) = s + , 
A = A - kizi - k2{l - zi) . 



Furthermore, the following relations hold 



-t 
-u 



Z2 



I-Z2 



A%M2 

Zi{l - Zi)' 

zi{s + fci') + (1 - z,)k2^ + 2A ■ fca, 

(1 - Zi){s + fci^) + Zik2^ - 2A ■ k2, 

[{l-z,)q-A]' + M^ 
[ziq + Af + M'' 

ZiV 



In this notation, the cross section ^2, M^, z/)for production of heavy quarks 

gluon collisions g*g* QQ equals (from pT|]) 

all{k,M,M\v) = -j^^yali^i^) j rf$(2)Z}(fci,fc2,z/, A,;2i,M2), 

and the matrix element D{ki, k2, ly, A, Zi, M^) is given by 



1 



{B + Cf 
2 '-,2 



iV, [(M2-n)(M2-t) ki'k 
1 / 1 



+ 



+ 



2 ((l-22)fcl' + (l-^l)fc2' + fci-fc, 



vs ki^k2^ 



where the following notation has been used 



^ ^ 1 {I - Zi){l - Z2)U Z/(l - 2:1 - 2:2) A-(fci + fc2) 

2 M2 - M 2s s ' 



1 ziZ2y z/(l - 2:1 - 2:2) A • (fci + ^2) 

2 ~ M2 -t 2i 1 



22 



References 

[1] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 71 (1993) 500; 71 (1993) 2396; 75 
(1995) 1451; B. Abbott et al. [DO Collaboration], Phys. Lett. B 487 (2000) 264. 

[2] J. Breitweg et al. [ZEUS Collaboration], Eur. Phys. J. C 6 (1999) 67; C. Adloff et al. [HI 
Collaboration], Nucl. Phys. B 545 (1999) 21. 

[3] J. Breitweg et al. [ZEUS Collaboration], Eur. Phys. J. C 12 (2000) 35. 

[4] J. J. Aubert et al. [European Muon Collaboration], Phys. Lett. B 106 (1981) 419. 

[5] C. Adloff et al. [HI Collaboration], Phys. Lett. B 467 (1999) 156; Erratum-ibid. B 518 
(2001) 331. 

[6] J. Breitweg et al. [ZEUS Collaboration], Eur. Phys. J. C 18 (2001) 625. 

[7] M. Acciarri et al. [L3 Collaboration], Phys. Lett. B 514 (2001) 19. 

M. Acciarri et al. [L3 Collaboration], Phys. Lett. B 503 (2001) 10; The OPAL Collabora- 
tion, OPAL Physics Note PN455. 

[9] S. Frixione and P. Nason, JHEP 0203 (2002) 053. 

[10] J. C. Collins and R. K. Ellis, Nucl. Phys. B 360 (1991) 3. 

[11] S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. B 366 (1991) 135. 

[12] E. M. Levin, M. G. Ryskin, Y. M. Shabelski and A. G. Shuvaev, Sov. J. Nucl. Phys. 53 
(1991) 657; 54 (1991) 867. 

[13] M. G. Ryskin, A. G. Shuvaev and Y. M. Shabelski, Phys. Atom. Nucl. 64 (2001) 1995. 

[14] A. V. Lipatov, V. A. Saleev and N. P. Zotov, Mod. Phys. Lett. A 15 (2000) 1727. 



[15] A. V. Lipatov, V. A. Saleev and N. P. Zotov, |hep-ph/0112TT4 
[16] H. Jung, Phys. Rev. D 65 (2002) 034015. 
[17] A. Szczurek, |hep-ph/0203050[ 



[18] C. B. Mariotto, M. B. Ducati and M. V. Machado, |hep-ph/0208"T55 



[19] G. Marchesini and B. R. Webber, Nucl. Phys. B 386 (1992) 215; S. Catani and F. Haut- 
mann, Nucl. Phys. B 427 (1994) 475. 

[20] B. Andersson et al. [Small x Collaboration], Eur. Phys. J. C 25 (2002) 77. 

[21] L. N. Lipatov, Sov. J. Nucl. Phys. 23 (1976) 338; E. A. Kuraev, L. N. Lipatov and 
V. S. Fadin, Sov. Phys. JETP 44, 443 (1976); ibid. 45 (1977) 199; I. I. Balitsky and 
L. N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 



23 



[22] M. Ciafaloni, Nucl. Phys. B 296 (1988) 49; S. Catani, F. Fiorani and G. Marchesini, Phys. 
Lett. B 234 (1990) 339; Nucl. Phys. B 336 (1990) 18; G. Marchesini, Nucl Phys. B 445 
(1995) 49. 

J. Kwiecihski, A. D. Martin and A. M. Stasto, Phys. Rev. D 56 (1997) 3991. 
H. Jung and G. P. Salam, Eur. Phys. J. C 19 (2001) 351. 

M. A. Kimber, J. Kwiecihski, A. D. Martin and A. M. Stasto, Phys. Rev. D 62 (2000) 
094006. 

M. A. Kimber, A. D. Martin and M. G. Ryskin, Phys. Rev. D 63 (2001) 114027. 
M. A. Kimber, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C 12 (2000) 655. 
K. Golec-Biernat and M. Wiisthoff, Phys. Rev. D 59 (1999) 014017; 60 (1999) 114023. 
N. Timneanu, presentation at the Small-x workshop, Lund, June 2002 
A. L Shoshi, F. D. Steffen, H. G. Dosch and H. J. Pirner, |hep-ph/ 0207287 . 
A. Gawron and J. Kwiecihski, [hep-ph/0207299 . 



N. Timneanu, J. Kwiecihski and L. Motyka, Eur. Phys. J. C 23 (2002) 513. 

V. Gribov and L. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438 and 675; L. Lipatov Sov. 
J. Nucl. Phys. 20 (1975) 94; G. Altarelh and G. Parisi, Nucl. Phys. B 126 (1977) 298; 
Y. Dokshitser, Sov. Phys. JETP 46 (1977)641. 

Y. L. Dokshitzer, D. Diakonov and S. I. Troian, Phys. Rept. 58 (1980) 269. 
M. Gliick, E. Reya and I. Schienbein, Phys. Rev. D 60 (1999) 054019. 
M. Gliick, E. Reya and I. Schienbein, Eur. Phys. J. C 10 (1999) 313. 
J. Kwiecihski, |hep-ph/0203T72 . 

M. Gliick, E. Reya and A. Vogt, Eur. Phys. J. C 5 (1998) 461. 
C. Adloff et al. [HI Collaboration], Phys. Lett. B 483 (2000) 36. 

V. M. Budnev, I. F. Ginzburg, G. V. Meledin and V. G. Serbo, Phys. Rept. 15 (1974) 181. 

R. D. Ball and R. K. Ellis, JHEP 0105 (2001) 053. 

T. Sjostrand et a/., Comput. Phys. Commun. 135 (2001) 238. 

H. Jung, Comput. Phys. Commun. 143 (2002) 100. 

S. P. Baranov, H. Jung, L. Jonsson, S. Padhi, N. P. Zotov, Eur. Phys. J. C 24 (2002) 425. 
M. Cacciari and P. Nason, Phys. Rev. Lett. 89 (2002) 122003. 

E. L. Berger, B. W. Harris, D. E. Kaplan, Z. SuUivan, T. M. Tait and C. E. Wagner, Phys. 
Rev. Lett. 86 (2001) 4231. 



24 



